I-balls 



S. Kasuya", M. Kawasaki'' and Fuminobu Takahashi^ 



" Helsinki Institute of Physics, P.O. Box 64, FIN-OOOI4, University of Helsinki, Finland 
Research Center for the Early Universe, University of Tokyo, Tokyo 113-0033, Japan 

(Dated: February 1, 2008) 

We find that there exists a soliton-like solution "I-ball" in theories of a real scalar field if the 
scalar potential satisfies appropriate conditions. Although the I-ball does not have any topological 
or global U{1) charges, its stability is ensured by the adiabatic invariance for the oscillating field. 
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I. INTRODUCTION 
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04 ' Scalar fields play important roles in theories of the early universe. It is believed that our universe experienced quasi- 
$H ■ exponential expansion phase (— inflation ) in its very early stage, which solves the flatness and horizon problems of the 
^ ■ standard cosmology and explains the origin of the density fluctuations of the universe such as observed by COBE jl| 

and other experiments H, y, Q . The inflationary universe scenario is realized by the vacuum energy of some scalar 
field (inflaton). After inflation, the inflaton starts to oscillate and decays into other particles which reheat the universe 
through thernialization processes. 

Similar dynamics is found in the Affleck-Dine mechanism for baryogenesis which is a promising candidate for 
^ ' explaining the matter-antimatter asymmetry of the universe. The mechanism makes use of a scalar fleld (AD field) 
00 corresponding to a flat direction in the scalar potential of the minimal supersymmetric standard model. During 
10 inflation the AD field has a large field value and oscillates when the effective mass becomes smaller than the Hubble 
CO ^ parameter after inflation. When the AD fleld starts oscillation, the baryon number is generated through the baryon 
i number violating term in the potential. 

' Recently, it was found that the oscillating AD field deforms into lumps of the scalar condensate called Q balls 

5 ■ SSI!- The Q baU is a non-topological soliton and its stability comes from the global charge (=baryon number) 
' conservation. The existence of the Q ball is crucial because it may significantly change the scenario of the Affleck-Dine 

baryogenesis fioj- The fragmentation into scalar lumps may also take place for the inflaton field. In fact, Enqvist et 
al [TJl pointed out that the oscillating inflaton fleld can fragment into Q balls. 
O 1' Since the Q ball is stable owing to the charge conservation, the scalar fleld responsible for the Q ball must be 
^ • complex. Then, a question arises whether or not a real scalar field deforms into lumps similar to the Q balls. At 
first glance, no stable lumps are formed because any conservation quantities like a global charge do not exist for 
the system of a real scalar field. However, the previous studies on the dynamics of scalar fields showed that some 
soliton-like objects are formed. For example, "oscillons" are formed for phase the double well potentials 0| and the 
axion field fragments into "axitons" 15]. In both cases the numerical simulations showed the existence of some scalar 
lumps inside which the scalar fields are rapidly oscillating. However, the reason why such quasi-stable scalar lumps 
are formed was not clear at all. Moreover, recently, McDonald 12] pointed out that in a hybrid inflation model the 
inflaton field can fragment into scalar lumps even if the scalar field has any conserved charges [see also Ref. [l3|]. 

Thus, it has been seen that real scalar fields fragment into quasi-stable lumps in numerical simulations for various 
situations. In this paper, therefore, we face the important problem concerning the real scalar dynamics, that is, 
what makes the scalar lump quasi-stable? Because the conserved baryon number plays a crucial role for stability 
of the Q ball, we need similar conservation quantity to stabilize the real scalar lump. In classical mechanics it is 
well-known that the adiabatic invariant exists for oscillating phenomena As will be seen later, we find that the 
adiabatic invariant can be extended to the field theories. Then, the existence of the stable lump can be explained by 
the adiabatic charge / ( see the next section for definition ) for the oscillating scalar field. We call this scalar lump 
"I-ball" , since the adiabatic charge / plays the same role as the global U{1) charges in the case of Q balls. We obtain 
the condition on the form of the scalar potential for the I-ball formation and derive the equation which determines 
the field configuration of the I-ball. In particular, it is found that the adibaticity requires the scalar potential to be 
dominated by a quadratic term. We also perform numerical simulations to confirm the existence of the I-balls for two 
types of simple potentials. 
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II. ADIABATIC CHARGE 



In this section we derive the conservation of an adiabatic charge, which guarantees the stabihty of the I-ball as 
discussed later. First we shortly review an adiabatic invariant in a mechanical system, according to Ref.|0|. Let us 
suppose that a system is executing a finite motion in one dimension and characterized by some parameter A(t) which 
specifies the properties of the external field. We assume that A varies slowly enough (i.e.,"adiabatically"): 



(1) 



where an overdot represents a derivative with respect to time, and T is the period of the motion. If A is constant, 
the system executes a strictly periodic motion with a constant energy E. For slowly varying A, the energy E varies 
slowly, while there is a quantity which remains constant, called an adiabatic invariant. This is written as 

hdira = ^<j>pdq, (2) 

where q and p are the coordinate and momentum, and the integral is taken over the variation of the coordinate during 
one period. Note that the periodicity plays an essential role in the proof of the existence for the adiabatic invariant. 
Before we go on to the case of the field theory, it will be useful to discuss a multi-dimensional system. Let us consider 
a system with any number of degrees of freedom {qi,Pi\, executing a finite motion in all the coordinates. It is assumed 
the variables are separable so that the action can be written as the sum of functions each depending on only one 
coordinate. As shown in Ref. 16|, the motion of the system is in general not strictly periodic, but conditionally 
periodic since the system passes arbitrarily close to a given state in the corse of a sufficient time. In fact, it is periodic 
only if the frequencies of all degrees of freedom are commensurable for arbitrary values of {qi}. In this case, there 
exists only one adiabatic invariant in the system. On the other hand, if the variables are not separable, there is no 
adiabatic invariant in general. However, if the Hamiltonian of the system differs only by small terms from one which 
allows separation of the variables, the properties of the motion are close to the periodic motion. 

Thus, in order to extend the adiabatic invariant to the case of the field theory, we have to impose the following 
two assumptions. First, the gradient energy is always sub-dominant everywhere, since the action of a scalar field is 
separable except the gradient term. Second, the scalar potential is quadratic, which would ensure the strictly periodic 
motion if it were not for the gradient term. In the following arguments, we adopt these two assumptions. 

Let us consider the system of a real scalar field ^, whose motion is finite and characterized by a parameter A(x, t). 
We suppose that A varies adiabatically, and that its dependence on the position x is weak enough as well. In the limit 
of constant A, the motion of the scalar field is homogeneous and periodic with the period T . The Lagrangian is given 
as. 



£ = -d^ct>d^<t>~V{ct>,\), (3) 



V{<i^,>)\t<^ = i^m'^\ (4) 

where we assumed that the potential is quadratic with the mass equal to m before the external field A(x, t) is turned 
on at t = 0. If A were constant, the energy-momentum conservation law would be given as 

d.T^'' - 0, (5) 

T^"' EE d^'(t)d''4)-r]^"'C. (6) 

We turn attention to its time component, 

d^.f = 0, (7) 

where an overdot represents 9t, and ?7^^ = diag(-f , — , — , — ). With slowly varying A for i > 0, the motion of the scalar 
field changes adiabatically. Then the energy-momentum current is no longer conserved: 

d,f^{d,\)^. (9) 



3 



If we average Eq. ^ over the period T, we have 



t+T 



dt 



'TJ, " d\ 

where the overhne represents the average over the period of the motion: 



djf" d(j) 



^4 



t+T 



dtZ. 



In deriving Eq. Hl()|l . we used the two assumptions, one of which states that A varies adiabaticaUy: 



(10) 



(11) 



(12) 



Therefore, A can be regarded as a constant, when § d(f> is integrated over the total motion of i 
fixed X. The other assumption is that A depends on the position weakly enough: 

VA 1 



A 



during one period at 



(13) 



where T can be interpreted as a typical spatial scale of the oscillating system. This condition is necessary because 
otherwise the gradient term becomes too large to be negligible. Note that the large gradient term also changes the 
value of 0^ significantly during one period, leading to the violation of the adiabaticity. For a small but nonzero VA, 
the motion is not strictly periodic, and the deviation from the orbit obtained as if A were constant can be estimated to 
be of the order of 54) ~ (f)TV\ after one period. We neglect such small corrections, and would like to focus attention 
on the leading terms in the following argument. Thus the motion of (p is approximated to be both periodic over T 
and homogeneous in the volume ^ = [x — T /2, x + T/2]. Hence there are four conserved quantities = jy T'^'^d^x 
along the path of the motion, and one can regard d^^cj) as a function of (0, J'^, A). 
By differentiating the equation : j'^{(j),dfj_(j>, X) = J'^ with A, we have 

dd„6 dX dX 



= 0, 



In deriving Eq. H14|) . we have used the relation 



dj^ _ 
'dX ^ ' 



^A~" 



(14) 



(15) 



where the non-diagonal components are higher order in A and VA, and hence we neglect them. With the use of Eq. 
Eq. l(Tn|) can be rewritten as 



a(9^</>) 

-1 



dX^^+^ 



9(9^0) 



(16) 



Hence we have 



d^i^j^dcbd^cfj ^0, (17) 

where the time component (/i = 0) is the leading term, and has a definite physical meaning. Thus, we are led to 
define the adiabatic charge / as 

^ d^x (18) 



2m 



where the factor 1/m in the definition here is introduced to make / dimensionless. Apparently / is conserved, 

dt 



(19) 



4 



III. CONDITIONS FOR EXISTENCE OF I-BALLS 



With the use of the adiabatic charge derived in the previous section, let us consider the condition that I-baUs are 
formed. We would like to focus attention on A in the first place. It specifies the properties of the external field, hence 
the energy of the system is no longer conserved for varying A (see Eq. 0). Alternatively, it could be argued that the 
role of A is played by the self-interaction of the scalar field, which must be such that the adiabatic conditions, H12(l 
and p3|l are satisfied. Then the total energy of the system including the interaction would be conserved, although 
the energy of the free part varies due to the self-coupling. In other words, there are two invariants, the energy and 
adiabatic charge for the whole system including 'the external field'. 

We assume that the scalar potential V{(j)) is given as 

V{4>)^^m^(t>^ + Vi{ct>), (20) 

where the self-interaction Vi{(j)) is small enough to respect the adiabatic conditions. That is to say, the adiabatic 
conditions are assumed to be satisfied by the dynamics of the system. First, we separate (j) into the rapidly oscillating 
part (P) and slowly varying part ($) as 

./)(x,t) =$(x,i)F(x,t). (21) 

For the general potentials, P might strongly depend on x, which leads to the large gradient energy (V0)^ and violates 
the adiabatic condition However, we can safely adopt the ansatz that P is homogeneous over a sufficiently 

large scale in which we are interested, since the periodicity of the system is guaranteed if the adiabatic conditions are 
satisfied. Thus, 

F(x,i)=P(0 (22) 

where P{t) oscillates between —1 and 1 with the period T. We also define 0(1) constants c„ ( n = 1, 2, 3, • • • ) for the 
later use. 

^ = c„$2", (23) 

02 ^ CiTO^^^ (24) 

where we used the fact that the scalar potential is dominated by the quadratic term. 

We take advantage of the method of Lagrange multipliers to look for the minimum of the energy E at fixed /, and 
minimize 

E^ = E + {d(i-^ f d^x4>^ 

= y^^xci Q|V$|^-^^$2-f y($)^ -H(w + m)/, (25) 



where w = w -t- m, and V{^) = V{(j))/ci. Assuming the bounce solution is spherically symmetric. 



E^ = j drA-Kr^ ci ^ 



1 ^ 



2\dr) ^m-^*') + (26) 



In order to minimize i?^ with respect to $, we have to seek for the bounce solution, which is equivalent to solving 
the equation, 

dr"^ r dr d$ ' 



where 



[/($) = r_$2_y(^)^ (28) 
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The bounce solution should satisfy the following boundary conditions [see Fig.^: 



dr 



0, 



■r=0 

$(oo) = 0. 

Hence, the bounce solution (I-ball solution) exists if the following inequalities are satisfied. 



$2 



< Lum < m 



(29) 
(30) 

(31) 



Now we investigate the constraint on the interaction Vi . For the field configuration to satisfy the adiabatic condition 
(|T!^ . we require \{d^/dr)/^\ <^ m. Since the spatial scale of the I-ball solution is determined by the mass scale of 
U{^), this requirement is equivalent to 



d2[/($) 




d^Vii'P) 


d $2 




d $2 



(32) 



where we used the inequality H31|l . If the interaction satisfies this constraint, the periodicity of the system is main- 
tained, and the I-ball configuration which minimizes the total energy of the system is also gently-sloping enough. 

Lastly, we comment on the cases that the scalar potential is not dominated by the quadratic term. Up to here we 
have shown that the adiabatic invariant can be found for a scalar field with the quadratic potential in the external 
field, and that the I-ball configuration minimizes the energy of the whole system including the interaction, which plays 
a role of the external field, for the fixed adiabatic invariant. Thus it is not evident from the preceding arguments 
whether the lumps are formed for other potentials. However, we performed the numerical calculations, and found 
that no quasi-stable lumps are formed for the several examples of such potentials. Hence this fact strongly suggests 
that the existence of I-balls is peculiar to the case where the quadratic term dominates the potential. 



IV. NUMERICAL SIMULATION 



Now that we have shown that the I-balls minimize the energy for fixed /, it must be then investigated whether 
such soliton-like objects are really formed. For this purpose we perform numerical calculation which follow the 
time evolution of the system. This is a quite nontrivial question, since the requirement that the system should 
vary adiabatically and its spatial distribution be gently-sloping enough ( "adiabaticity condition" ) must be satisfied 
dynamically. As concrete examples we take both the gravity-mediation like potential and m^^^ _ ^4 potential. 

First let us suppose the potential is given as 

V{cb) = Im'f (^l + Xlog(^^)) , (33) 

where is a renormalization point to define the mass, and the K term is the one-loop correction which is assumed 
to be negative. The I-ball equation Eq. l|77|l reads 

^2$ 2d$ , , / 

^ + -^+-o^'J'--^'J>^log(^j=0, (34) 

where Wq is defined as 

ujl = Ljm-m^{l + K). (35) 

From numerical calculations, it is suggested that a Gaussian ansatz is a reasonable approximation to the I-ball solution 
for this potential. If we insert the Gaussian ansatz, 

$(r) = $(0)e-'''/^' , (36) 

into the I-ball equation Eq. (|34(l . we obtain 

4r2 6 \ / , 2i^m2 „\ ^ ^ 
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where we have set = $(0)/2. Thus we see that the same form is obtained in the first and second terms. This 
requires that 

Ri - (38) 



which roughly corresponds to the inverse of the most amphfied mode fcmax, and this correspondence is checked 
numericaUy [see Fig. |2]. Also the maximum growth rate of the linear perturbation is estimated to be of the order 
\K\m in the similar way as in the case of Q balls 0, 0|. Hence the adiabaticity condition is satisfied for the 
gravity-mediation like potential. 

We perform the numerical simulation to confirm whether the I-balls are formed with the radius obtained above. 
For the numerical calculation, we take the variables to be dimensionless as follows. 

<P = — , 
m 

r mt, 
Xi ^ mXi, 

Hi = kjm, (39) 
where ki is a wave number in the Xi direction. The initial conditions are taken as 

ip{0) = 1.6 X 10^ +(5i, 

V'{0) = -I + S2, (40) 

where the prime denotes the derivative with r, and we set K = —0.1. (5's are fluctuations which originate from the 
quantum fluctuations, and their amplitudes are taken to be 10^^ times smaller than the homogeneous mode. We have 
confirmed that the smaller fluctuations just delay the formation of I-balls. First we check that the adiabatic condition 
is satisfied up to the linear growth of perturbations. Fig. |21 shows the numerical result of the instability band for the 
initial conditions stated above. It can be seen that the instability band roughly corresponds to the inverse of the 
radius Ri obtained analytically. We present the result of numerical simulation in two dimensional lattices in Fig. |31 
from which one can see that the energy density p deforms into lumps, identified as I-balls. Also Fig. 0] represents 
the profile of the scalar field <I> inside the I-ball, and the analytic solution obtained above agrees quite well with the 
numerical result, which suggests the Gaussian ansatz is appropriate. Though there is some deviation in the outer 
region, it is irrelevant since the absolute value of the scalar field is much smaller than that at the center of the I-ball. 



Next we take the following potential. 



n^) = ^-v - + (41) 



where we have added the ip^ term for the stability of the vacuum. Here a is positive to satisfy the condition H31|l . 
Actually, for negative a, no I-balls are produced in the numerical simulations. As in the previous case, we take the 
variables to be dimensionless as Eq. (|39|l . The initial conditions are taken as 

<p{0) ^ 1.0 + 5i, 

^'{0) = O + S2. (42) 

We set a = 0.1, b — 0.005, and &' amplitudes are taken to be 10^^ times smaller than the homogeneous mode. As 
long as —(j)'^ term is much smaller than the mass term, the growth rate a and the instability mode k are given by 

««.(|^)*N. 

3 
2' 



< k2 < |a$2 . (43) 



Some reasons can be adduced; (1) The adiabatic conditions might not be satisfied well around the surface, since the homogeneous mode 
is relatively small compared to the fluctuations. (2) Actually the I-balls are not isolated. 
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Since we take = m as an initial condition, a(kmax) and ^/kl^ are much smaller than m, so the adiabaticity 
condition is satisfied up to linear perturbation. The numerical result of the instability band for the initial conditions 
stated above is shown in Fig. 13 We can see that the instability band coincides exactly with that obtained analytically. 

We present the result of numerical calculation in two dimensional lattices in Fig. |^ The energy density p deforms 
into I-balls in this case, too. We have performed numerical simulations with different a, and found that all of these 
results generally look alike. When the I-balls are newborn, they are almost spherically symmetric, and the gradient 
energy is subdominant everywhere, suggesting that our formulation is valid. As the system evolves in time, their 
shapes deviate from the spherically symmetric one, and become irregular. Finally, they decay into random phases 
where the kinetic, potential and gradient energies are all same order. The lifetime of the I-balls r becomes shorter 
and they are formed earlier as a increases, but its typical value is O{10^m~^) for a = 0.1. This is consistent with the 
facts pointed out in Refs. [T^IT5l |. that the oscillons and axitons have very long but finite lifetime. Hence, it is certain 
that the I-balls for this type of potential decay in the end. The decay might be induced by the small deviation from 
the adiabaticity of the dynamics, which induces the decoherence of the oscillating scalar field inside I-balls. But it 
needs further investigations to make clear how the decay proceeds. 

V. CONCLUSION 

We have studied the system of a real scalar field and found the solution of the quasi-stable scalar lump, I-ball. The 
stability of the I-ball can be explained by the adiabatic invariant charge /, which does stem from the dynamics of 
the system, not any symmetries. For the I-ball solution to exist, the scalar potential should be dominated by the 
quadratic term {nn?(tP') and satisfy the condition Ij^lf) which is almost same as that for the Q-ball. 

Furthermore, we have performed numerical simulations and have found that the quasi-stable I-balls are really 
produced and their properties are in agreement with the theoretical predictions. Since scalar fields prevail in theories 
of the early universe, the I-balls may be formed and play important roles in various cosmological processes . 
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FIG. 3: Spatial distribution of the energy density p for the gravity-mediation like potential. Inside the I-balls, the scalar field 
cj> oscillates rapidly. 
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FIG. 6: Spatial distribution of the energy density p for the m^fj)^ — (f>* potential at time corresponding to r = 450, when the 
I-balls are newborn. Inside the I-balls, the scalar field (p oscillates rapidly. 



